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! Abstract 



Phenomenological studies support the applicability of naturalness and naive 
dimensional analysis to hadronic effective lagrangians for nuclei. However, 
one-baryon-loop vacuum contributions in renormalizable models give rise to 
unnatural coefficients, which indicates that the quantum vacuum is not de- 
scribed adequately. The effective lagrangian framework accommodates a more 
general characterization of vacuum contributions without reference to a Dirac 
sea of nucleons. 
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Variations of the Walecka model are widely applied in descriptions of nuclear structure 
and reactions llj-Q. These models are appealing both for their phenomenological success 
and for the simplicity and economy of the physics at the mean-field level. Most practical 
applications used for quantitative nuclear structure phenomenology include contributions 
from valence nucleons only; the framework is often labeled a "no sea" approximation to 
indicate that the Dirac sea of nucleons is neglected p|. 

The approach to the relativistic nuclear many-body problem known as quantum hadro- 
dynamics (QHD) was originally based on renormalizable field theory so that systematic 
calculations with a finite number of parameters are possible, at least in principle. In partic- 
ular, one-baryon-loop vacuum effects can be included as a way to add vacuum dynamics to 
the "no sea" physics. These effects have a natural interpretation as the response of a filled 
Dirac sea of nucleons to the presence of valence nucleons The resulting "relativistic 
Hartree approximation" (RHA), however, does not provide an acceptable description of the 
properties of finite nuclei, at least by the standards of modern successful mean-field models 

The deficiencies of the RHA could imply that one simply needs to work harder to describe 
vacuum dynamics within the QHD framework. It is known, however, that including higher- 
order corrections in a simple loop expansion makes the phenomenology worse ^ , so that any 
improvement (if it exists) would require significantly more complicated diagrams, such as 
those involving vertex corrections or short-range correlations |TO| , p!I| . Moreover, one would 
like to structure the hadronic many-body framework so that the appealing, intuitive, and 
successful mean-field theory is indeed a good starting point. This suggests an effective field 
theory (EFT) approach based on the "modern" viewpoint of nonrenormalizable lagrangians 
||T^^T7[. This viewpoint is well known in the particle physics community (for systems at 
zero density) and is adopted in applications of chiral perturbation theory to meson-meson, 
meson-nucleon, and nucleon-nucleon scattering. However, it is not well known or applied 
in the community of nuclear physicists who use relativistic models for nuclear structure. 

From the EFT perspective, a description including only valence nucleons and classical 
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("mean") meson fields can still incorporate vacuum dynamics, hadron compositeness, and 
many-body correlations, albeit approximately [p|, p^pO| . The price to be paid is that all 
possible terms consistent with the underlying symmetries (excluding redundancies) must 
appear in the effective lagrangian. Nevertheless, by relying on the concept of "naturalness" 
(as defined below), it is possible to systematically truncate the effective lagrangian, leaving 
only a finite number of unknown parameters; moreover, recent fits to empirical nuclear 
properties using this framework give strong evidence that the model parameters are indeed 



natural ai8-20 



From this point of view, the RHA in renormalizable QHD models is simply one specific 
prescription for determining an infinite number of parameters in the effective theory, namely, 
the coefficients in the scalar effective potential and of terms involving derivatives of the bo- 
son fields. Here we assess the relevance of the RHA prescription by examining the size of 
these coefficients. We find that the RHA leads to unnaturally large coefficients, in disagree- 
ment with results obtained from fits to empirical nuclear properties This implies that 
although it may be possible to explicitly include vacuum dynamics by calculating baryon 
vacuum loops, it is much more efficient to include them implicitly in the small number of 
natural parameters contained in the truncated effective lagrangian. 

There have been more formal criticisms of the RHA: (1) the RHA vacuum contributions 



violate Nc counting rules motivated from quantum chromodynamics (QCD) pl| , |22| , (2) 
the RHA neglects the compositeness of the nucleons p3|- p5|] , and (3) the treatment of the 
vacuum in terms of A^A^ pairs alone is simply wrong, or at best incomplete [p^ , ^5| -|28|l . It 



seems much more compelling to us that one should avoid explicit baryon-loop calculations 
simply because the empirical properties of nuclei show that dynamical vacuum effects are 
quite modest and can be described with a few adjustable parameters. 

Georgi and Manohar p9| , p!5| have proposed a naive dimensional analysis (NDA) for as- 



signing a coefficient of the appropriate size to any term in an effective lagrangian for the 
strong interaction. This NDA has been extended to effective hadronic lagrangians for nu- 
clei |l^j2^. The basic assumption of "naturalness" is that once the appropriate dimen- 



sional scales have been extracted using NDA, the remaining overall dimensionless coefficients 
should all be of order unity. For the strong interaction, there are two relevant scales: the 
pion-decay constant ^ 93 MeV and a larger scale 0.5 ^ A < 1 GeV, which characterizes 
the mass scale of physics beyond Goldstone bosons. The NDA rules prescribe how these 
scales should appear in a given term in the lagrangian density: 

1. Include a factor of I/Ztt for each strongly interacting field. 

2. Assign an overall normalization factor of /^A^. 

3. Multiply by factors of 1 /A to achieve dimension (mass)^. 

4. Include appropriate counting factors (such as l/n\ for 0"). 

The appropriate mass for A might be the nucleon mass M or a non-Goldstone boson mass; 
the difference is not important for numerical assessments of naturalness, but will be relevant 
for the N(. counting arguments considered later. 

As an example of the NDA prescription, a term in the scalar effective potential takes the 
form 

1 / A\n 



The coupling constant k„ is dimensionless and of 0(1) if naturalness holds. Until one can 
derive the effective lagrangian from QCD, the naturalness assumption must be checked by 



fitting to empirical nuclear data. Such fits give strong support for naturalness |p!8| , p0 |. 

We can assess the naturalness of the RHA vacuum contributions by matching the results 
in a renormalizable model to an effective mean-field theory in which the Dirac sea degrees of 
freedom are excluded by construction. The matching to the RHA is conveniently made with 
an effective action, in which the Dirac sea contribution is easily isolated. (See, for example, 
Ref. ip.) The one-loop Dirac sea effective action V in models with Yukawa couplings to the 
nucleon takes the (unrenormalized) form 

r' = -th Tr ln(^7^9^ - M* - g^^V^) , (2) 



where M* = M — gs(f), 4> and are neutral scalar and vector fields, and are their 
couplings to the nucleon, and the trace is over spatial and internal variables. For time- 
independent background meson fields, the effective action is proportional to the energy.[] 

The expression in Eq. can be renormalized and evaluated as a derivative expansion 
[|§ in the and V^" fields: 

This expansion in inverse powers of M* converges rapidly for finite nuclei [pIH 3^. The 
effective potential Ucs and the coefficient functions Zi can be further expanded in (infinite) 
polynomials in (p with well-behaved coefficients; that is, there is a local expansion of F', 
which can be absorbed into an effective lagrangian for nuclei. Thus the finite Dirac sea 
contribution can be reproduced by an effective lagrangian treated at the mean-field level 
with only valence nucleons, as long as all possible terms (generally nonrenormalizable) are 
included. 

The effective potential Ucs{4') is found by evaluating the trace in Eq. (|^) with constant 
fields. This expression is divergent and must be regularized and renormalized to obtain a 
finite result, which takes the general form [Q] 

^eff(0) = --^(M*Mn — + ^a„M^-"(^?,0)") , (4) 

where 7 is the spin-isospin degeneracy and the a„ are dimensionless constants. The scale 
parameter fi is typically chosen to be M. The counterterms ao and ai are fixed by requiring 
Ugs to be zero and a minimum in the vacuum (0 = 0). The others are fixed by prescrip- 

^The contribution to the energy from Eq. (^) can be written more transparently as a sum over 
single-particle energies of occupied states in the filled Dirac sea Q. In nuclear matter, this (un- 
renormalized) energy density is AS = V~'^ Ek,A[(k^ + M*2)V2 _ (^2 _^ jvf2^i/2j_ 
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tion. Note that in a renormalizable model, only the first four powers of are available as 
counterterms. In an effective lagrangian, however, all powers are present. 

The most common prescription has been to choose the a„ to cancel the first four powers 
of (p appearing in the expansion of the logarithm |Q. One finds in this case (for 7 = 4) 

1 r d) 7 13 9^ n 

U,sU) = — t- (M - $)4 ln(l - — ) + - -M^^^ + — M<l>3 - — $4 
<=ttvr; 4^2 [V ' ^ 2 3 12 J 

_M^r$5 $6 <l>^ 4!(n-5)!$" ^ 

where $ = (7s0- When (p (or M*) is determined at each density by minimization, Ucs{4') is 
the finite shift in the baryon zero-point energy that occurs at finite density and is analogous 
to the "Casimir energy" that arises in quantum electrodynamics. 

To evaluate the size of the one-loop vacuum correction, we apply the NDA. Based on the 
scahng rules discussed above, a term of 0(0^) should be scaled as 

where we have associated A with the nucleon mass M. (See, however, the comments on 
Nc scaling below.) If this contribution is natural, any residual overall constant should be of 
order unity. However, if we perform a similar scaling on the leading term in Eq. (|]), we find 

4ir= 5AP 6 /J ^ \5\fy J ' ^ ' 

where we used 47r/^ ^ M and Qs ~ M/f^^. Thus the one-baryon-loop contribution to the 
vacuum energy is roughly two orders of magnitude larger than naturalness requires. It is 
not hard to show from Eq. (^ that all higher powers of $ contain essentially the same large 
overall factor. 

We can make an instructive comparison of natural and unnatural coefficients using the 
linear sigma model, generalized to include a neutral vector meson. Variations of this model 
have been used to investigate the role of chiral symmetry in nuclear structure |3^|1|,^ .f\ After 



^The limitations of tfiis approach are discussed in Refs. |^ and l9|. 
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following conventional procedures to introduce a nonzero expectation value for the scalar 
field and then shifting the field, we obtain the model lagrangian (with = for simplicity): 

l^au, = i'ihi.d'' - g^i^V^" - (M - g^a) - ig^'j^r • Tr]^/' 

+ ^{d^ad'^a - mla^) - ^{d^V, - d.V^f + ]^mlV^V^ 

+ \d,n . d^n + g.^a{a^ + - gl^,{a' + , (8) 

where the associations a ^ (p, g.,, ^ gs, and m^j — > nis should be made for our discussion. 
Comparing to Eq. (|l|), with A identified as ms, we find K3 = — K4 = —3, so that the nonlinear 
parameters are natural at the mean-field level.0 However, if one includes the one-baryon- 
loop vacuum corrections, renormalized in a fashion that preserves the chiral symmetry []35|,p|, 
one finds unnatural corrections to the cubic and quartic couplings: Ak,^ = 2M^/7r^/^ ^ 20, 
= — SM^/vr^/^ ^ —80. The quintic and higher corrections are exactly the same as in 
Eq. (D. 

These unnatural coefficients generate correspondingly large corrections to the conven- 
tional linear Walecka model mean-field theory (MET). If we adjust the model parameters 
to reproduce the standard nuclear matter properties in the RHA (equilibrium at Eermi mo- 
mentum /Cp = 1.30 fm^^ with a binding energy of 15.75 MeV), the baryon effective mass at 
equilibrium becomes M*/M 0.73. This translates into a change in the scalar potential 
$ from 430 MeV in the MET to 250 MeV in the RHA. This is a large effect, particularly 
since it is the (p^ term that makes the difference. When a general effective lagrangian is fit 
to nuclear properties, terms of this order play essentially no role ||20[ . 



Contributions to the scalar potential at equilibrium density from Ues{4>) and from the 
original tree-level lagrangians are plotted in Eig. |l| for the Walecka model and in Eig. |^ for the 
linear sigma model. In each figure, the crosses indicate contributions to the energy density 
estimated using naturalness, with error bars allowing for a range of A and M* /M. A steady 



^This is to be expected; the model builds in at tree level a realization of QCD chiral symmetry 
and Goldstone-boson physics, which is the same physics behind the NDA. 
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FIG. 1. Contributions to the scalar potential per particle in nuclear matter from the n^'^-order 
terms of the form <I>" for the RHA model. The crosses are estimates based on Eq. The arrow 
indicates the total binding energy eo = 15.75 MeV. 

decrease with increasing n is evident, which motivates the truncation of effective lagrangians 
with natural coefficients. M* /M in natural models (and phenomenologically successful ones) 
is typically between 0.60 and 0.66. For those values, the RHA 0($^) contribution would 
be as large as a typical 0($'^) contribution in a natural model (see Fig. |I|) and would 
prevent a successful description of nuclear structure. In fact, the unnatural and unbalanced 
0($^) term drives M* /M to its self-consistent value of 0.73. Higher-order contributions are 
essentially negligible, because of the natural factor of (p/f-^ (and a combinatoric factor) that 
accompanies each higher power. In the linear sigma model, which has unnatural 0{^"^) and 
0($'^) contributions from the RHA, M*/M is driven to 0.9. 

It is possible to devise a prescription |^ that leads to an effective mass M* /M and a 
nuclear compressibility consistent with a reasonable (although not optimal) fit to properties 
of finite nuclei (see Ref. []19[ for the criteria). However, this requires choosing a„ coefficients 
in Eq. (^) to achieve sensitive cancellations between terms of different order in the effective 
potential, so as to neutralize the effect of the unnatural (p^ contribution. The unnaturalness 
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FIG. 2. Contributions to the scalar potential per particle in nuclear matter from the n^'^-order 
terms of the form for the linear sigma model plus RHA. The crosses are estimates based on 
Eq. (|^). The arrow indicates the total binding energy eo = 15.75 MeV. 

of the vacuum nucleon loop contributions is generally characterized by the changes in the 
coefficients that accompany 0(1) changes in the scale /i [say, from /i to e/i in Eq. (Q)], 
as in Ref. [|29|. The subsequent changes in the a„ coefficients are obtained by expanding 
7M*^/(47r^). These changes are large compared to the natural size implied by Eq. (|l]). 

Based on the strong empirical evidence for naturalness, we conclude that the treatment 
of the quantum vacuum at the one-baryon-loop level is, at best, inadequate. Although the 
concept of a Dirac sea is compelling for nuclear physicists because of the analogy to the Fermi 
sea, the explicit calculation of these effects prejudices the description of the vacuum dynam- 
ics and (to date) has not yielded results consistent with nuclear structure phenomenology. 
Moreover, the self-consistent, valence- nucleon-only theory is covariant, causal, and internally 
consistent by itself, and the empirical evidence shows that the vacuum degrees of freedom 
can be included implicitly by a small number of local interactions among the mesons and 
valence nucleons. Note that the omission of explicit dynamical contributions from the Dirac 
sea does not mean that one can discard negative-energy solutions entirely; they must be 
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retained to ensure the completeness of the Dirac wave functions in certain calculations of 
density-dependent effects (for example, linear response) [0,^ . 

It should also be emphasized that an effective lagrangian allows for a more general 
characterization of the vacuum dynamics than that arising from baryon loops. The explicit 
calculation of counterterms in the effective lagrangian is unnecessary, since the end result is 
simply an infinite polynomial in the scalar field, with finite, unknown, and apparently natural 
coefficients arising from the underlying dynamics of QCD. To have predictive power, one 
must rely either on the truncation scheme provided by naturalness (see Ref. @), so that 
only a small, finite number of unknown coefficients are relevant, or on some other dynamics 
to constrain the form of the renormalized scalar potential. For example, a simple model is 
used in to show how the broken scale invariance of QCD leads to dynamical constraints 
on the scalar potential, and fits to the properties of finite nuclei also generate coefficients 
that are natural. 

Another potential difficulty in the explicit calculation of baryon vacuum loops is that 
RHA vacuum contributions violate counting rules motivated by consistency with QCD 
| 2T| , |22[] . The expected Nc scaling of the coupling for an n-meson vertex is O^N^^"^^"^) 



The Nc scaling property of a vertex in an effective lagrangian or the RHA is established 
with the associations M oc 0{Nc), rris cx 0(1), and gl oc 0{N^ |^. [Note also that 
/tt oc 0{NI/'^)] Thus one may simply inspect the coefficient of 0". For the RHA, the result 
is 0{N^~"'^'^), which exceeds the expected scaling by a factor of A^"^ pT] , p2 . 



In contrast, the effective lagrangian approach, with our definition of naive dimensional 
analysis and naturalness, is consistent with the scaling expected from QCD if one associates 
A in the 0" terms with nis (and not M). For example, applying the Nc scaling rules to the 
coefficient of 0^ produces 

-5|/.^A^| -5^(m3V^)| ^ 0(1) X 0(Ar-3/2) ^ 0(^-3/2) 

which is consistent with the usual Nc counting. It is important here that the normalization 
factor A^/^ comes from the meson mass term, with the meson mass rris being 0(1) rather 
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than 0{Nc). Thus Nc counting imphes that the vacuum response is qq in nature; the response 
of a A^A^ vacuum is unhkely to agree with the large Nc hmit. 

Finally, we comment on the statement often made or implied in the literature that 
hadronic field theories treat nucleons as point particles and so do not account for the com- 
positeness of nucleons. The use of effective lagrangians with fields for composite particles 
is by now well established (e.g., chiral perturbation theory with nucleons and pions). In 
Ref. [^, the electromagnetic structure of the nucleon manifestly appears at the mean- field 
level using a hadronic effective lagrangian for nuclei. The key features are the inclusion of 
nonrenormalizable interaction terms and the use of a derivative expansion to incorporate 
the nucleon compositeness; indeed, the allowed freedom in the definition of the boson fields 
shows that even renormalizable terms like 0^ and 0^ implicitly include the effects of nucleon 
substructure. Thus, the deficiencies of the RHA are not intrinsic to the use of nucleons (and 
the Dirac equation) to describe nuclei, but arise because the implied vacuum dynamics is 
incorrect or incomplete. 

In summary, we have examined the relativistic Hartree approximation (RHA) to renor- 
malizable hadronic field theories in the context of modern effective field theory. The param- 
eters obtained in phenomenologically successful mean-field models of finite nuclei exhibit 
naturalness, as anticipated from naive dimensional analysis of the strong interaction. How- 
ever, the parameters implied from vacuum loops in the RHA are unnatural. Since the 
phenomenological parameters in successful models implicitly include the effects of vacuum 
dynamics, we conclude that the explicit treatment of the vacuum in the RHA, which involves 
the Dirac sea of nucleons, is inadequate. In contrast, a natural effective model with only 
valence nucleons, including all possible (nonredundant) terms, is consistent with nucleon 
compositeness, Nc counting, and nuclear structure phenomenology. Extensions to include 
higher-order vacuum loops will be discussed in subsequent work. 
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